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State estimation is usually analyzed in the situation when copies are in a product state, either 
mixed or pure. We investigate here the concept of state estimation on correlated copies. We 
analyze state estimation on correlated A'' qubit states, which are permutationally invariant. Using a 
correlated state we try to estimate as good as possible the direction of the Bloch vector of a single 
particle reduced density matrix. We derive the optimal fidelity for all permutation invariant states. 
We find the optimal state, which yields the highest estimation fidelity among the states with the 
same reduced density matrix. Interestingly this state is not a product state. We also point out 
that states produced by optimal universal cloning machines are the worst from the point of view of 
estimating the reduced density matrix. 
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I. INTRODUCTION 

In the standard approach to quantum state estimation, 
one is given a state which is a product state of N iden- 
tical copies of a quantum system in an unknown quan- 
tum state. The task is to find the optimal measurements 
on the N copies and the optimal guessing strategies in 
order to learn what is the actual state of the quantum 
system. Optimal strategies differ depending on many fac- 
tors: What figure of merit describes the accuracy of the 
guess? What is the dimensionality of the Hilbert space 
of a single copy? How many copies of a state are given? 
What is the prior knowledge about the state (in partic- 
ular is it mixed or pure)? Is there a limitation on the 
allowed measurements (for example restriction to local 
measurements - separate measurement on each copy)? 

Many results have been obtained in recent years on the 
subject. In particular the optimal fidelity for estimating 
the state of N qubits in an unknown pure state was found 
in 0. More generally, optimal fidelity for estimating 
arbitrary pure state of a d-level system was derived in 
0, 111 • The case of mixed states is more demanding, and 
the general results for the optimal estimation fidelity are 
known only for the case of qubits nil. Estimation via 
local measurements was analyzed in the case of qubits 
for pure and mixed states H, 0] . 

In all the references mentioned above, state estimation 
is performed on the input state of the form: p®^, where 
p is either pure or mixed state of a single copy. This is 
a product state in which copies are not correlated. The 
intriguing questions is: what is the optimal state estima- 
tion when copies are correlated? Do correlated copies al- 
low for a better estimation of the reduced density matrix, 
or on the contrary, do correlations diminish the estima- 
tion fidelity? In other words is the product state p^^ 
optimal from the point of view of state estimation? 

In this paper we analyze the state estimation on corre- 
lated qubits. We analyze correlated density matrices of 
N qubits, which are permutationally invariant, and have 
single particle density matrices p, described by a Bloch 



vector of length r. We pose the problem of the optimal 
estimation of the direction of the Bloch vector given the 
correlated state of N qubits. We find out that the prod- 
uct state is not the optimal state for the purpose of 
state estimation. We derive the optimal state, and the 
corresponding optimal fidelity of estimation of the opti- 
mal state. Additionally we also point out the fact that 
the state coming out from the optimal cloning machine 
is in fact the worst state from the point of view of state 
estimation, among all the states with the same single 
particle reduced density matrix. 

II. CLONING AND STATE ESTIMATION - THE 
MOTIVATION 

There is an interesting connection between optimal 
cloning and optimal state estimation, which appeared to 
be fruitful in deriving optimal fidelity for N — > Al cloning 
of qubits, using the known result for state estimation 
Interestingly, benefits flew also in the opposite direction, 
as knowledge of the optimal local cloning fidelity of a 
pure state of a d-level system(9| allowed the derivation of 
the optimal estimation fidelity in this case 0] ■ The main 
result relating cloning and estimation is, that the local 
fidelity of — > oo optimal cloning of a quantum state is 
equal to the fidelity of optimal estimating the state using 
N copies. 

In the derivation of the above result, as one of the 
steps in the proof, the state estimation on the clones 
coming out of the cloning machine is considered. Clones 
coming out from the optimal N ^ M cloning machine 
are always correlated (with the exception of the trivial 
case N = M). This is an example of the situation where 
one has to consider state estimation on correlated copies. 
For the sake of the proof authors make use of the fact 
that the state of the clones is supported on the symmetric 
subspace, and thus can be decomposed as Q 

p = ^a,|^.)(V.|®^, (2.1) 
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where a; = 1 and at can be both positive or negative 
(from now on p denotes the fuU density matrix, while 
p denotes single particle reduced density matrix). This 
decomposition allows one to relate the fidelity of state es- 
timation on correlated clones to the fidelity of estimation 
on uncorrclatcd pure states 

Some additional observations can be made concerning 
the relation between cloning and state estimation. Con- 
sider optimal universal N ^ M cloning of qubits. If the 
input A'^ copies were in the state \'ip)®^ then the single- 
particle reduced density matrix of each clone is given by 

Ilia: 

=K^,M)|^)(^| + i^^l^^l, (2.2) 



where 



r{N, M) ^ 



N M + 2 
M N + 2' 



(2.3) 



When the number of clones M goes to oo the reduced 
density matrix of each clone reads: 



N 



Pn ^ , 



7V + 2 



1 



-1. 



(2.4) 



Estimation of the initial state \ip) can be done using these 
infinite number of clones (in reality a very large number), 
and the resulting fidelity is f3| : 



F ■ 



N + 1 
N + 2' 



(2.5) 



This result is obtained using the fact that the full density 
matrix of clones is supported on symmetric subspace, and 
this fidelity is just the fidelity of optimal state estimation 
on N qubits. 

Let us assume now (counter-factually) that after N 
oo cloning, the reduced density matrix of each clone is 
given by the same formula as in Eq. (|2.4() . but unlike in 
real cloning machine, the full density matrix of clones is 
a product state: pprod — (Pn ^ 00)*^°°- What is the fi- 
delity of estimating the initial state \tp) using the state 
Pprod- Each single-particle density matrix, has a Bloch 
vector pointing in the direction lip). Since we have in- 
finite number of such density matrices independent one 
from another, with arbitrary high accuracy we can esti- 
mate the direction of the Bloch vector, and consequently 
the state \ip). The fidelity of estimating \ip) using the 
state Pprod is thus 1. Consequently, if clones were really 
in such a product state we could beat the state estima- 
tion fidelity limit. This observation shows, that the re- 
strictions on cloning coming from the state estimation are 
not confined to the restriction on local fidelity of cloning. 
State estimation imposes also some necessary degree of 
correlation between the clones. 

In this context, it is interesting, to study in general the 
state estimation fidelity on correlated copies. One could 
then impose limitations on possible states of clones solely 
using limits from state estimation theory. 



The other problem, on which the analysis of estima- 
tion on correlated copies could shed some light, is re- 
trieving information from quantum systems after they 
have passed through a noisy channel. Consider for ex- 
ample N qubits in a pure state lip) that pass through a 
noisy channel, and each qubit at the output has reduced 
density matrix p, with the Bloch vector shrunk by cer- 
tain factor. One can ask the question, what kind of noise 
leads to smallest (largest) loss of information about ini- 
tial state \ip)7 Is the noise acting independently on qubits 
less harmful than the one inducing correlation between 
them? 



III. STATE ESTIMATION ON CORRELATED 
COPIES 

The problem of state estimation on correlated copies 
can be stated as follows. Consider a density matrix p 
describing the state of N d-level systems. In order not 
to distinguish any of the systems we shall assume that p 
is permutation invariant i.e.: 



npnt = p, 



(3.1) 



where 11 is an arbitrary permutation of N systems. Sin- 
gle particle reduced density matrices of each system are 
equal and denoted hy p = 1^2. ..nP- The goal is to find op- 
timal measurements on the state p and optimal guessing 
strategies in order to estimate the state p. The measure- 
ment is described by a positive operator values measure 
[llj |: -P/i = 1- With each measurement result a guess 
Pf_i is associated. If F{p^,p) denotes a figure of merit, 
describing how close a guess is to the actual single 
particle state p, then the average fidelity of the state es- 
timation is given by the formula: 



i^ = ^F(p,„p)Tr(P^p), 



(3.2) 



where averaging is performed over the ensemble of un- 
known initial states p. 



IV. OPTIMAL STATE ESTIMATION ON A 
CORRELATED N QUBIT STATE 

We now move on to the specific case of state estimation 
on a correlated state of N qubits. 

Single particle density matrix of each copy is given by: 



P = r\xp)(ip\ 



-1 



1(. 



(4.1) 



We assume that the length r of the Bloch vector r is fixed 
and our goal is to estimate its direction (this is for exam- 
ple the case when relating cloning with state estimation). 
In other words we want to estimate the pure state lip), 
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given a mixed state p with single particle reduced den- 
sity matrices p. In what follows we shall assume that we 
have no prior knowledge about the direction of the single 
particle Bloch vector. Averaging in Eq. H3.2I) will thus be 
performed over the ensemble of input states of the form; 

~p{U) = [/«^pt/t®^ (4.2) 



, where U are elements of the SU(2) defining representa- 
tion, distributed uniformly with respect to the invariant 
measure on the SU(2) group. 

As a figure of merit we shall adopt the fidelity: 
F(/9^, |^)(V'|) = (V'IP/ilV')- If the full density matrix is 
a product one: ppiod = then the optimal fidelity of 
state estimation given in can be rewritten as: 



= E Xl^^J. + i) + ^f'^' H3 + 4,) -!] + (!- rr-^] , (4.3) 

j=s 



where dj (for j ^ N/2) is defined as: 



dj = 



N 



N 



(4.4) 



1, and s = (1/2) for N even (odd). 



U'N/2 

Instead of the product state, consider now the state 
Psym, which is supported on the symmetric subspace (for 
example this is the case of the clones coming out from 
the optimal cloning machine), and has the same single 
particle reduced density matrix p. The two states Pprod, 
Psym coincide only when r = 1 (pure state case). The 
optimal fidelity of estimating the direction of the Bloch 
vector f (|4.1|l reads: 



F -1(1 



N 



N 



(4.5) 



This result reveals what is known from analyzing cloni_nj 
in stages, namely that shrinking factors multiply 
Note that here we abstract from the problem of cloning, 
and analyze general state Psym supported on symmetric 
subspace, which may be not be obtainable by means of 
optimal cloning machines. The proof of the above for- 
mula, using the methods from goes as follows. The 
fidelity of estimation can be written as: 



F = (V'IE;.PA.Tr(P^p)||V') 



(4.6) 



(V' I Post I V"), 



where post is the state we estimate on average (average 
with respect to different measurement results fi). The 
state estimation procedure may thus be described by a 
completely positive map £: 



Pest ^£{p), 



(4.7) 



Thanks to the choice of the ensemble of states over which 
we average, the search for the optimal state estimation 
procedures can be restricted to the set of covariant oper- 
ations £ i.e.: 



£{U^^pU^^^) =U£ip) Ul 



(4.8) 



Covariance property of £ implies that the action of £ on 
a pure product state yields: 



#>(0I 



1 - a 



-a, 



(4.9) 



where a > and does not depend on |0) (this is because 
the result must commute with any rotation around the 
axis defined by the Bloch vector of 4>). 

Using the decomposition (|2.1|l of the matrix p, and the 
fact that 



P 



= Tr2...w(p)=Ea^|V'»)(V'»|, (4.10) 



the fidelity of state estimation can be written as (averag- 
ing can be dropped out thanks to covariance): 



F 



= (^1 



1 



-1 



1^) 



2 

(4.11) 

The highest fidelity corresponds to the highest possible 
value of a. From Eq. H4.9|l one sees that the highest pos- 
sible a corresponds to the estimation strategy optimal for 
pure product states i.e. a = N/{N + 2) l\. Substituting 
this value into the Eq. H4.11|) one gets (|4.5|) . ■ 

Interestingly the relation -Fprod ^ -fsym holds, with 
equality only when r = or r = 1. This may not seem 
obvious from the formulas (|4.3|) . H4.5|) . but will become 
clear in Sec. Ivl 

We now derive the optimal estimation fidelity for a 
more general class of density matrices, namely when p is 
an arbitrary permutation invariant density matrix. As a 
byproduct of the derivation we will be able to write down 
the state Popt, which allows for the best estimation of the 
direction of the single particle Bloch vector, among all 
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the density matrices p with fixed single particle reduced 
density matrix p. Interestingly this in general will not be 
the product state Pprod- 

We shall follow the methods presented in |5|. In the 
space of N qubits (spin 1/2) one can introduce an or- 
thonormal basis |j, m, a) [j — s, . . . , N/2 denotes the 
total spin, whin s = (1/2) for N even (odd)] such 
that: (i) for a fixed j and a vectors \j,m,a) span the 
(2j + 1) dimensional space, in which an irreducible rep- 
resentation of SU(2) acts, (ii) for a fixed j and m vectors 
|j, TO, a) form a basis for dj dimensional H4.4|l irreducible 
representation of the symmetric group (Sn)- This repre- 
sentation correspond to a Young diagram with two rows 
of lengths N — 2j and 2j. Let Uj^a be unitary matri- 
ces that connect spaces with fixed {j, a) with the space 
(j, 1): \j, m, 1) = Uj,a\j, TO, a). For convenience it is cho- 
sen that: |j, m, 1) — \j, to) ® |vI'_)®^/2^J, where \ j, m) is 
the symmetric state of 2j qubits with the total spin j and 
the projection of the spin onto the r direction (direction 
of the single particle Bloch vector) equal to, while l^*-) 
is the singlet state. 

The above decomposition of the Hibert space into sub- 
spaces invariant under permutations and SU(2)'^^ is a 
manifestation of the Schur-Weyl theorem ■ The sym- 
metric group acts irreducibly in the subspaces with fixed 
j and TO. This allows us to write down the most general 
permutationally invariant density matrix p, in the simi- 
lar way that the product state p*^^ was written down in 
i: 

N/2 ^ d, 

P^Y.Pi^Y.Pj^'-^ (4.12) 

j — s a—1 

where 



where 

. _ ^ 3 + m R _ J - ^ y 

m——j m——j 

The Bloch vector corresponding to the density matrix pj 
reads: 

r- = [C, + C*,i{C, - C*),A, - B,] (4.15) 
and its length is equal to: 

r, = \r,\ = ^[A,-B,Y-2\C,\^ (4.16) 

Using formulas (|4.12|) . 14.13|l . we obtain the single par- 
ticle reduced density matrix p: 

P = Tr2,...,^(p) = Y.P^ f^p, + ^^l) (4.17) 

Comparing it with the single particle reduced density ma- 
trix H4.1|l . one arrives at the following relation between 
local Bloch vector f and vectors fj : 

2 

This condition guarantees that the state p, parameterized 
hy Pi, and Xi,m,m', yields the single qubit reduced density 
matrix with the Bloch vector equal to f. 



P],a - Uj^aPj.lUj^^ 

= p,^\^^){^^r^/^~l ^^^^^ 
3 

Pj = E ^],m,m'\j,m'){j,m\. 
m^m' — —j 

Coefficients pj are nonnegative and satisfy Yli^llPj — 
while for a fixed j the matrix Xj^m,m' is a positive semi- 
definite matrix with trace one. Summation over a assures 
that the density matrix p will indeed be permutation in- 
variant. 

Let us now express the single particle reduced density 
matrix p with the help of coefficients Pj, Xj^m.m' ■ First 
we calculate the single particle reduced density matrix of 
the symmetric matrix pj. Written in the basis 
it reads: 



Pj = Tr2,...,2, (p,) = ) (4.14) 



A. Optimal estimation strategy 

We are now prepared to derive the fidelity of the op- 
timal state estimation performed on the state p parame- 
terized hypi, Xj^m,m' ■ First let us observe that similarly 
as in 0, ^3 1 W6 do not loose optimality if we start with 
projecting the state p on a subspace with given (j, a). 
This is simply a measurement projecting on subspaces 
with respect to which p is already block diagonal, and 
therefore this operations cannot decrease accessible in- 
formation on the state. If the measurement yields a ^ 1 
then Uj,a is applied in order to move the state to the sub- 
space (j, 1), and consequently move singlets to the last 
N/2 — i qubits. Since singlets yield no information on the 
direction of the local Bloch vector they can be discarded. 

After this preprocessing, with probability pj, we have 
at our disposal a 2j qubit state pj. Further procedure 
depends on our knowledge on the structure of correlations 
in the state p. 

Knovirn structure of correlations. Knowing the 
structure of correlations means that we know everything 
about the state p except for its possible local rotation 
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with {7®^ . This means that even though we do not know 
the direction of the local Bloch vector f, we know the 
relative orientation of vectors r} in Eq. H4.18(l . In other 
words, we know all scalar products between r and Vj. 

This information can be helpful in increasing our fi- 
delity of estimation. The state pj which we obtain from 
the preprocessing stage is a 2j qubit state supported on 
the fully symmetric space. We will use the estimation of 
the direction of its local Bloch vector Vj (optimal esti- 
mation of symmetric state described at the beginning of 
this section) as an information about the direction r. We 
may ask if it is advantageous to make additional trans- 
formation of our guesses about the direction fj in order 
to make them closer to the vector r (higher overlap) and 
hence increase the fidelity. Note that we do not know 
the direction r so this transformation has to be defined 
independently of r. 



We know the scalar product ' 



which tells us that if 



the direction rj is set to point to the north pole then flics 
somewhere at a parallel of colatitude 9j = arccos(fj • r). 
Since we completely do not know at which point on the 
parallel lies the vector f, then thanks to the rotational 
symmetry of the problem, it may shown that if 9j < 90° 
any transformation of guesses of fj will actually lead to a 
lower average overlap - it is optimal to keep the guesses 
as if we wanted to guess vector fj at the north pole. If, 
however, 9j > 90° then it is clearly preferable to make 
a reflection with respect to the center of the sphere as if 
we wanted to estimate the vector pointing to the south 
pole. 

Consequently, the optimal way to use the state pj for 
determining the direction f, is to estimate the direction 
of its local Bloch vector fj and in the case 9j > 90° make 
an additional reflection of the guesses [20j |. 

We can introduce p'j which is equal to the reduced 
density matrix pj of the state p^, whenever 9 < 90°, and 
when 9 > 90° its Bloch vector is additionally reflected. 
The optimal fidelity of estimating the direction of f using 
the state pj reads (compare the derivation l|4.11ll 'l: 



2j 



2J , 



1 



2j 



a 1^) 



max(^j , Bj 



2j + 2' 



(4.19) 



where \4>) is a pure state with the normalized Bloch vec- 
tor f/r and max{Aj, Bj) appears as a result of the con- 
ditional reflection of the Bloch vector of p'j . 

Expressing Fj explicitely using Bloch vectors, and 
summing over j with weights pj we get the fidelity of 
estimation when the structure of correlations is known: 



-^known 2 



N/2 



(4.20) 



Unknown structure of correlations. In this case 
we do not know anything about the structure of correla- 
tions in the state p. In particular we do not know rela- 
tive orientations of vectors fj, and cannot judge whether 



9j < 90° or 9j > 90°. As a result we do not know when 
to apply reflection, and the natural |2l| approach in this 
case is not to apply any reflection at all. Therefore, we 
estimate the direction of fj using the state pj , and regard 
this result as an estimate of the direction of f. 

The only difference in the formula for the fidelity as 
compared to Eq. (|4.20|) is the lack of the absolute value, 
since we do not apply reflection if 9j > 90°. The fidelity 
reads: 



N/2 



F 



unknown 



^.7 + 1 



(4.21) 



3=s 



Obviously, in general we will arrive at a lower estimation 
fidelity than the one given in Eq. (|4.20fl . 



iV QUBIT STATE OPTIMAL FOR STATE 
ESTIMATION 



We are now able to answer the question posed at the 
beginning of this paper. With a fixed length of the Bloch 
vector r of a single particle reduced density matrix, what 
is the state p optimal for estimation of the direction of 
the Bloch vector? Mathematically the goal is to opti- 
mize estimation fidelity Fknown or i^unknown under the 
constraints: 



N 



N/2 

E 



3=s 
N/2 

^Pj^l, Pj>0, rj<l 



(5.1) 



(5.2) 



Since -Fknown or -Funknown depend on fj only via scalar 
products fj ■ f , we can always replace vectors fj , by their 
projections on the direction f. Both the constraint 1)5. 1|) 
and the fidelities are not affected by this. From now on, 
we can thus regard all vectors fj lying in the same line 
as f, and the constraint H5.1|l reads: 



N/2 



N 



(5.3) 



where ± corresponds to vector fj being parallel and an- 
tiparallel to f respectively. 

Let us start with the case when structure of correla- 
tions is known. 

Knovirn structure of correlations. Looking at 
Eq. H4.20|l . we notice that it is always advantageous for 
the value of fidelity to take maximal possible length of 
vectors fj, i.e. rj = 1. We can always do it, and still 
keep the constraint (|5.3|) satisfied, since if necessary we 
can choose some vectors fj to be antiparallel, and change 
weights Pj appropriately. Consequently, we have to max- 
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imize the quantity (note jr*, ■ r|/r = 1): 



^ known 



N/2 



(5.4) 



Notice that j /{j + 1) is an increasing function for posi- 
tive j. We should therfore try to give as much weight as 
possible to pj with the highest j. Ideally we would like 
to have Pn/2 = 1 all other pj — 0. This, however, is 
only possible when r = 1. In other cases the condition 
(|5.1|) would be violated. If r < 1 we need to include some 
other Pj , the optimal choice is to take additionally only 
the second highest, i.e. with j — N/2 — 1. In order to sat- 
isfy the condition (|5.1|) , and maximize A the direction of 
rjv/2) should be respectively parallel and antipar- 

allel to f (in terms of Xj^m,m' coefficients, this means that 

Xn/2.N/2,N/2 — li Xn/2-1. -{N/2-1). -{N/2-1) — ^: ^ud aU 

other coefficients are zero). 
Constraint H5.1|l implies that 



- {PN/2N/2 ~ {N/2 - 1)) = r. (5.5) 

Together with the fact that pn/2 +PN/2-1 = 1 we get 
explicitely: 



PN/2 
PN/2-1 



N{r + l)-2 

2N -2 

N ~rN 

2N -2 



(antiparallel) 



and the corresponding estimation fidelity reads: 



7V2 



{N -l){N + 2)' 



(5.6) 



The fidelity depends very weakly on r, which means that 
even if we have very mixed local density matrices - very 
short local Bloch vector r - we can have a very good esti- 
mation of the direction r. This can be easily understood, 
when realizing that the short Bloch vector f arises from 
subtraction of two long vectors (see Eq. (|5.5|) '). Moreover, 
we can perfectly distinguish subspaces with j — N/2 and 
j — N/2 — 1. Hence, when we measure j = N/2 — 1 we 
know that 7^/2-1 points in the opposite direction than r, 
and we refiect it. As a result we get very good estimation 
of the vector r. In the asymptotic limit ^ cx), up to 
the leading order in 1/7V the fidelity reads: F = 1 — 1/iV. 
It does not depend on r and is the same as the asymptotic 
fidelity when estimating N copies of a pure state. 

For this strategy to be applicable, we need to perform 
collective measurements. If we restrict ourselves to lo- 
cal measurements assisted with classical communication, 
we can not discriminate between spaces with different j 
with high confidence. This leads to an interesting situa- 
tion when global measurements allow for almost perfect 
estimation of direction, while local measurement, even in 
a limit of large number of copies yield very poor estima- 
tion quality |l4| . 



Unknown structure of correlations. Our goal is 
to find a state that under the fixed length r of the local 
Bloch vector will maximize the fidelity in Eq. H4.21|l . As 

a matter of fact we need to maximize: 



^unknown 




(5.7) 



where ± indicates whether the vector r*, is parallel or 
antiparallel to r. The presence of the ± sign in Eq. 1)5. 7|) 
in contrast to Eq. (|5.4|l . makes a great difference when it 
comes to finding the optimal state for state estimation. 

Notice that because of the constraint H5.3|l . the sum of 
enumerators in Eq. I|5.7(l is fixed. For a moment, let us 
consider only the situation when all terms in Eq. (|5.7|l are 
positive. The fact that the denominator increases with j, 
together with the fixed sum of enumerators makes it not 
desirable to have contribution from terms with large j. 
On the other hand, for high values of r it impossible fulfill 
the constraint (|5.3|) without such terms. As a result, 
intermediate j, which are as low as possible, yet high 
enough to fulfill (|5.3|) . will prove to be optimal. 

Observe, that it is not desirable to have rj < 1, at 
least for j > s. If for any j' > s we have rj> < 1, then 
it is better to take r-,' = 1, decrease pj' , and increase pj" 
for some j" < j'. This will increase Eq. H5.7|l and at the 
same time allows to keep the constraint (|5.3f) satisfied. 

Let us proof the following lemma, still assuming that 
terms in Eq. H5.7|l are positive: 

Lemma. If for a given r the optimal fidelity correspond 
to pj> = 1 for certain j' > s, then after decreasing r by 
a small amount e (0 < e < 2/A^), the optimal fidelity is 
achieved by taking the parameters pj of the form: pjr — 
l — eN/2, Pj'-i = i-N /2, and the remaining equal zero. 

Proof. From Eq. (|5.1|l we have j' — Nr/2. First we 
allow only one additional pj" [j" ^ j') to have a nonzero 
value. After decreasing r by e, constraints are the fol- 
lowing: pj> + Pj" = 1, Pj'j' + Pj"j" = N{r — e)/2 and 
j" < j'. The Aunknown now reads: 



Ne 



^unknown 



J 



2(/-/')V/' + l J' + I 



.] 



f + 1 
(5.8) 

The above formula attains the highest value for the high- 
est possible j" i.e. j" — j' — I. Allowing only one ad- 
ditional Pj we find that the optimal case corresponds to 
j" = j' — 1. Thanks to the simple (linear in pj) form of 
Eqs. H5.7I5.3|I . fidelity for any other combination of pj, 
with more than two nonzero terms, can be written as 
a convex combination of fidelities corresponding to case 
with only two non zero pj. Consequently, the optimal 
case is that with only two nonzero terms i.e. pji , Pj'-i. 
Using the constraint 1)5. 3|l we arrive at: pj' = 1 — eA^/2, 
Pj'-i = eN/2. The discussion is only valid for e < 2/N . 
For e > 2/N , pj' according to the derived formulas would 
be negative. ■ 

If r = 1 , there is no freedom for pj H5.3|) , and we have to 
choose Pisi/2 = 1. This choice corresponds to the product 
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of pure states p = When we lower r by a 

small amount r = l — e(0<e< 2/iV), then using the 
lemma, the optimal state is given by: Piv/2-1 = -^(1 — 
r)/2, pt,,2 = 1 - N{1 - r)/2. For r = 1 - 2/N we obtain 
the situation in which only Pn/^-i = 1 is nonzero, and we 
can again use the lemma. We can carry on this process 
and find optimal states for smaller r till we reach a state 
where ps = \. 

It remains to discuss, whether negative terms in 
Eqs. H5.7I5.3|I can increase the fidelity. At first sight, 
it seems that they cannot. Taking negative term, means 
that in order to fulfill H5.3|l , we have to use positive terms 
with higher j , which will effectively make A smaller. This 
reasoning is not valid, however, for highly mixed local 
state (very low r). 

Consider first the case when N is odd. If r = f/A^, 
then according to the above discussion, the optimal state 
has pi/2 = f. For r < 1/N we cannot use the above 
scheme, since j — 1/2 is the smallest j. One of the 
options is to shorten ri/2 in order to fulfill (|5.3|) . but a 
better strategy is to use a negative term with j = A^/2, 
and keep ri/2 = 1- Because of j + 1 in the denominator 
of Eq. (|5.7() . we gain more by keeping ri/2 = 1, than we 
lose due to negative contribution from j = N/2 term. 
This also explains why we use j — N/2 and not a term 
with some other j - by doing so we subtract as little 
as possible, as the denominator of j = N/2 term is the 
largest. 

For N even, consider r = 2/N, for which the optimal 
state corresponds to pi = 1. Lowering r, we can keep the 
constraint (|5.3|) by using j = term. Although this is the 
optimal strategy according to the lemma, it is no longer 
optimal if we allow negative terms. This is because, the 
term with j = is the only one which does not contribute 
to the fidelity at all. It is more advantageous to keep pi as 
large as possible, and fulfill the constraint H5.3|) by adding 
negative contribution from j — N/2, then we gain more 
thanks to a larger pi term than we lose from subtracting 
j ~ N/2 term (again thanks to the denominator j + 1). 

Summarizing this lengthy reasoning, the optimal A^- 
qubit state for state estimation with fixed local Bloch vec- 
tor length r, when structure of correlations is not known, 
has all Vj = 1 and coefficients pj read: 

• For r > S/N. 



rN 



Pj-i = J 
For r < S/N. 

PS/2 

Pn/2 



1 

rN 

'Y 



when j 



rN 

1 < < 7 

- 2 -•' 



{r + l)N 

N + S 
S-rN 

N + S 



(antiparallel) 



where S = 2 for A^ even, and S" = 1 for A odd. 

With the above expressions one can easily calculate 
the optimal state and using Eq. H5.7|) the corresponding 
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FIG. 1: The figure presents the comparison of the optimal 
fidelity of state estimation on four different 6-qubit states: the 
state supported on the symmetric subspace (solid curve) , the 
product state p®^ (dashed curve), the state which we have 
found out is optimal from the point of view of state estimation 
when the structure of correlations is not known (dotted curve) 
and the optimal state when the structure of correlations is 
known (dot-dash curve) The goal of estimation is to estimate 
the direction of the Bloch vector of single particle reduced 
density matrices. The value of fidelity is plotted against the 
length of the Bloch vector r. For fixed r the state supported 
on the symmetric subspace (the kind of state that comes out 
from the optimal cloning machine) yields the lowest value of 
fidelity. 



optimal fidelity of estimation of the optimal state. In 
the asymptotic limit A^ — > 00, up to the leading order 
in \/N the fidelity reads: F = \ — l/{rN), and depends 
explicitely on r in constrast to the case when correlation 
structure is known. 

In Fig. n we present the plots for the optimal estima- 
tion fidelity on the state of six qubits in four cases; when 
p is: Psym (supported on the symmetric subspace), Ppiod 
(product state), Punknown (the optimal state for estima- 
tion when the structure of correlations is not known), 
Pknown (the Optimal state for estimation when the struc- 
ture of correlations is known). The length of the Bloch 
vector of the reduced single particle density matrix r is 
depicted on the x axis. The edgy shape of the dotted 
line is due to the fact, that in the case of unknown struc- 
ture of correlations, different pairs of pj are nonzero for 
different values of r. 

The 

Psym state appears to be the worst among these 
four states, from the point of view of state estimation. 
Actually, it is the worst state in general for estimation, 
when the structure of correlation is known. This can be 
seen from the Eq. (|4.2U|) . Only Pm/2 = 1 is non zero in 
this case, so only the term with the biggest denominator 
contribute. This makes the fidelity the smallest possible 
among all state with the same r. This in particular proves 
the inequality i^prod > -Fsym mentioned in Sec. IIVI 

Notice also that -Fknown, and -Funknown do not converge 
to 1/2 as r ^ 0, contrary to Fgym and Fpiod- Even if local 
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state are very mixed (arbitrary close to maximally mixed 
state), the estimation fidelity can be significantly higher 
than that 1/2 (which corresponds to random guessing). 

The state supported on the symmetric subspace psym, 
is interesting as this is the kind of state coming out from 
the optimal cloning machines. This is no coincidence that 
the state coming out from the optimal cloning machine 
is the worst from the point of view of state estimation. 
When we are doing M — > iV cloning, we can produce 
at the output only such a state of N copies that can 
be estimated with fidelity no higher than the optimal 
state estimation on M pure state. Otherwise, we would 
violate state estimation bounds. In cloning we insist on 
the highest possible fidelity of clones. We attempt to 
obtain the highest possible value of r (which determines 
local fidelity of clones). Doing so, and in order not to 
violate state estimation limits, we are forced to end up 
with the state which is the worst from the point of view 
of state estimation, among all the states with given r. 

The problem we have considered here, carries some 
similarities to the problem of the optimal encoding and 
decoding of a spin direction 

El EE III El m. In the 

optimal encoding and decoding problem, one considers 
N spins, and tries to find the optimal states of N spins 
that would be optimal for the purpose of transmitting 
the information about a direction in space. The optimal 
state for this task is in general {N > 2 El|) not a prod- 
uct states of N spins, but a correlated state. This shows 
that the information can be better encoded in correlated 
state than in a product one. Notice, however, that this 
result could not have been directly translated into our 
approach, as in our case we have strong restrictions on 
allowed states. We have considered permutationally in- 
variant states, and the optimal state (class of states) for 
state estimation that we have found, are states optimized 
under the constraint of fixed single particle density ma- 
trices. 

We should also stress, that we have considered the 
problem of the estimation of the direction of the Bloch 
vector of the single-particle density and not the single- 
particle density matrix itself. In other words, we have not 
attempted to estimate the length of the Bloch vector, but 
rather assumed it is known. Hence, there remain an open 
problem of the optimal estimation on correlated copies 
of the single-particle density matrix, when the length of 



the Bloch vector is unknown. This approach would be 
much more demanding, and in particular it is not clear 
what kind of distribution of correlated density matrices 
one should consider. We leave this problem open, yet we 
conjecture that also in this approach one should arrive at 
the conclusion that the optimal state for state estimation 
would be a correlated rather than a product state. 



VI. SUMMARY 

In this paper we have analyzed the problem of state es- 
timation done on correlated copies in order to determine 
the direction of the Bloch vector of the single-particle 
density matrix. In the case of qubits we have found the 
optimal fidelity of estimation on A'^-qubit states which 
are permutationally invariant. The state which is the 
optimal from the point of view of state estimation was 
derived. The optimal state is in general not a product 
state but has some correlations between qubits. As 
a result we have found out that the correlated state can 
be a better source of information about the single particle 
reduced density matrices than a product state. In some 
cases, though, large correlations prevent us from extract- 
ing information about single particle reduced density ma- 
trices, as it is the case with clones coming out from the 
optimal cloning machine. One is not entitled, however, 
to claim in general that correlation always worsen our 
ability to extract information about single particle den- 
sity matrix, as the optimal state for state estimation is 
in fact correlated. 
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